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Abstract: We consider the 2D quenched-disordered q— state Potts ferromagnets and show that at self- 
dual points any amalgamation of q — 1 species will fail to percolate despite an overall (high) density of 
1 — Further, in the dilute bond version of these systems, if the system is just above threshold, then 
throughout the low temperature phase there is percolation of a single species despite a correspondingly 
small density. Finally, we demonstrate both phenomena in a single model by considering a "perturbation" 
of the dilute model that has a self— dual point. We also demonstrate that these phenomena occur, by a 
similar mechanism, in a simple coloring model invented by O. Haggstrom. 

1 Introduction 

The purpose of this note is to address the question of whether there are "natural" translation and 
(lattice) rotation invariant ergodic measures on configurations rj £ {0,1} Z , d > 2, for which site 
percolation occurs when the density of occupied sites, p — (rji), is very close to zero and/or fails to 
occur if the density is close to one. It is known that neither phenomenon occurs for systems with 
product measures (independent percolation) [151 117] . where, in most situations, there is a sharp 
percolation transition when p^ 1 is of the order of the lattice coordination number (e.g. 2c?). The 
question then is what happens when the occupancy of different sites is specified in some natural way 
which is not independent. Of particular interest are equilibrium spin-systems and non-equilibrium 
stationary states of interacting particle systems such as the voter model, the contact process and 
their various generalizations 23J. Furthermore, one may also consider projections of such measures, 



1 



i.e. given a measure on [f2] z , where Oi £ ft is the set of possible states at i £ Z d , we define an 
occupation variable rji = 0, 1 , with rji = 1 if G Cl\ and rji = if <Xj £ fi\fii. Well known examples 
are the fuzzy q-state Potts model with q = r + s and rji = 1 if 0$ G (1, ..,r), and systems where 
a.; el, and 7jj = 1 if > /i and zero otherwise [TTJ [TO1 [5], 

Adams and Lyons pQ considered these types of question for measures on homogeneous trees 
that are invariant under the graph automorphisms. They found that whenever the density is large 
enough there indeed is percolation, see also [18] and [7]. Such a result does not hold for Z d ; it is 
easy to construct counterexamples which exploit, in an obvious way, the vanishing surface to volume 
ratio of regular sets [15]. Also, for d = 2 one may consider the infinite cluster of supercritical 
independent percolation near threshold and simply declare the complimentary sites to have rj\ = 1; 
these complimentary sites will have density close to one and fail to percolate. Higher dimensional 
analogs may even be possible, albeit at different thresholds (c.f. [2]) and the converse phenomenon 
can also be constructed by these means: Consider a slightly super-critical percolating cluster and 
declare its compliment to be vacant (rji = 0). While such examples are not particularly interesting 
in their own right, they demonstrate the existence of mechanisms for these phenomena which might 
occur in more realistic systems. 

Along these lines, a more interesting set of examples concern the so-called divide and color (DaC) 
models invented by Haggstrom [20]. As we will show, for these models in d= 2 there are no density 
limitations for the absence or presence of percolation. But even the DaC models are a bit artificial 
and one should ask whether there are more natural examples where such phenomena can occur. 

In this note, the principal system under study is the 2D disordered Potts ferromagnet. First 
we will demonstrate that at self-dual critical points, there is an absence of percolation despite a 
density which may be arbitrarily close to one. The complimentary result, namely percolation at 
small density would follow immediately from continuity of the magnetization - a currently open and 
challenging problem. (We will state some inconclusive results concerning this question.) However, 
almost trivially, the diluted versions of the disordered Potts models exhibit this property throughout 
the low temperature phase (i.e. when the magnetization is positive) provided the media itself is near 
threshold; we will provide a formal proof. Finally we will combine the two sets of results via the 
consideration of models with strong and weak interactions. These models inherit enough of the 
features of the dilute model to display the small density percolation property yet, unlike the dilute 
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models, can have a self-dual point. 

Most of the above stated results are proved in the context of a fixed realization of the disorder. 
However, since any particular realization is not translation invariant we do not satisfy the naturalness 
criterion that was stated at the beginning of this section. Therefore we will consider the entire 
problem from the perspective of the quenched measures. These objects are manifestly translation 
invariant and we will demonstrate various other amenable properties, such as strong mixing. 

The remainder of this paper is organized as follows: In section 2 we define the DaC models 
and, thereafter, we state and prove our results for this system. In section 3, we define the general 
quenched disordered Potts ferromagnets but with the emphasis on d — 2. We define the quenched 
measures for these systems; certain general properties of these quenched measures are stated (with 
the proofs postponed). We then prove our main results: (i) At certain self-dual critical points, 
any one of the species is sufficient to prevent percolation of the combined efforts of all the others, 
(ii) For the dilute model, with the active bonds in slight excess of the percolation threshold, there 
is a non-trivial low temperature phase characterized by the percolation of just one of the species 
notwithstanding (for g> 1) the paucity of the overall density of this species. Finally we construct 
a weak-bond/strong-bond disordered model with a self-dual point and, at least in the vicinity of 
the transition temperature, a regime of percolation at low density. In the next section, we will 
provide some discussion and finally, in the appendix, we prove the stated properties of the quenched 
measures in the disordered Potts systems. For further background discussion, see [15j Section 5.2. 

2 Critical 2D DaC Models 
2.1 Definitions 

On the "graph" Z 2 let p 6 [0, 1], q an integer > 2 and let Ai, . . . X q be in (0, 1) with ^ X q = 1. The 
relevant measure, which we denote by [i p .\ (with A = (Ai, . . . A g )) is defined as follows: Each edge 
is occupied (or vacant) independently with probability p (and (1 — p)). For each given configuration 
of edges, the finite connected clusters of sites are independently "colored" in one of q ways with 
respective probabilities Ai, . . . X q . If the bonds do not percolate, then the model has been defined. If 
there is percolation, p > p c = i, one can define q distinct "extremal" measures, /ip'(Ai...A g ), where 
the infinite cluster is assigned the jth color and convex combinations may be chosen according to 
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any desired prescription. The extreme measures can be justified/constructed via the limit of finite 
volume arrangements wherein the boundary of the region is deemed to be of the particular color 
chosen. 

Despite its somewhat artificial appearance, the DaC model has some interesting features (not 
discussed in the present work) and in fact can actually be realized in a number of circumstances. In 
particular, for Ai = • ■ • = X q it is the zero temperature limit of the random diluted g-state Potts 
model that will be featured in the next section. Further, if one constructs the g-state diluted voter 
model, with an initial distribution that is a priori equal and independent it is not hard to see (e.g. by 
the methods of [23] Chapter V) that the infinite time limit of this setup is precisely the DaC with 
corresponding p and Ai = • • ■ = X q . 

2.2 Statements and Proofs 

Our first theorem, the key parts of which are actually contained in 20J, reads as follows: 

Theorem 2.1 For all Ai £ (0,1), the measures ^ Pc ,(x 1 ,x 2 ) have, with probability one, no (connected) 
site percolation or *- connected site percolation for either of the species. 

Remark 1. We remark that the case of interest in the present work is when one of the A's is small; 
the setup in |20j was the completely symmetric p = X\ = A2 = ^. We further remark that by the 
obvious method of aggregating colors, the above extends to versions of the DaC with more colors. 
In particular, with q species, the union of any of the q — 1 species is denied percolation by the single 
remaining color. 

Proof. The critical bond percolation model on Z 2 has the key feature that there are infinitely many 
closed circuits of occupied bonds surrounding the origin and infinitely many closed circuits of dual 
bonds - edges on (Z + i) 2 - traversal to the vacant bonds surrounding the origin [17]. Since the 
latter separate the former, it is evident that there are infinitely many disjoint and separated occupied 
bond circuits surrounding the origin each of which belong to clusters of color type 1 or color type 2 
with respective probabilities Ai and A2. This rules out the possibility of the origin belonging to an 
infinite connected or even * -connected cluster of either type. □ 

*-percolation means an infinite cluster under the relaxed rules that neighbors or next nearest neighbors are 
considered connected. 
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The above theorem is complimented by the situation above threshold: 

Theorem 2.2 For every e > there is a region of p > p c and Ai such that Mwa! a 2 ) ^ as P erc °lation 
of the first color but the overall density of the first color is less than e. 

Proof. Consider the measure Mp(A t a 2 ) w ^ n P ( on ly slightly) in excess of p c . It is clear, by condi- 
tioning on the events that a particular site does and does not belong to the infinite bond cluster, 
that the color 1 density is given by Poo(p) +Ai(l — Pooip))- Here Poo(p) denotes the fraction of sites 
belonging to the infinite cluster of the bonds. Since, for the 2d bond problem, P oc (p) is known to 
vanish continuously [25] the result follows. □ 



3 Disordered Potts Ferromagnets 
3.1 Background Considerations 

We let {Ji.j | i,j neighbors in Z 2 } denote a collection of iid non-negative random variables with 
common distribution F not concentrated at a single point; we emphasize that we may allow Jjj = 
with non-zero probability, a case we will refer to as dilution. The bond-random Potts model is 
defined, for fixed realization of the (Jij), by the formal Hamiltonian 

where denotes nearest neighbors, J represents the collection (Jij) and <7i € {1, . . .q} are the 
spin-variables associated with the sites of Z 2 . Most often, we will absorb the usual temperature 
parameter into the definition of J. 

For fixed J in a finite volume A, the Gibbs measures with specified boundary conditions are 
defined in the usual fashion. E.g. if A C Z 2 and the boundary OA has fixed boundary spins ag^, we 
may write 

p — HjCcaIcsa) 

Pj?a>a)= - 7 * 3A (3-2) 

where H$(a\ | o-qa) denotes the object in Eq. (|3.1[) restricted to i and j in A U <9A, a\ and ag\ are 
notations for the collections of spins on these respective sites and the spins o~q\ are to be considered 
as fixed, Zj°£ is the partition function which normalizes P. Of exclusive interest in this work are the 
one ([!]) and free ([f]) boundary conditions. The former case is defined by simply setting all of oq\ 
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to be in the first spin state. (Since the Hamiltonian is completely symmetric under permutations of 
spin-states, all results concerning [l]-boundary conditions apply equally well to the (q — 1) others 
that are defined similarly.) The latter is described by the stipulation that Jjj is set to zero whenever 
ieA and j G <9A. Infinite volume thermodynamics is defined via the partition function in the usual 
way. All the standard thermodynamic functions (and their derivatives with respect to external 
fields) emerge as almost sure quantities which are independent of boundary conditions; c.f. [16] and 
references therein. 

For finite volume, the quenched measures are defined by averaging the probabilities in Eq. (|3.2[) 
over all realizations of the interactions according to Yin j) dF(Jij) (written as dF) - where the 
product runs over the appropriate set, usually i and j in A U dA). These averaged objects will 
be referred to as the finite volume quenched measures and denoted by Q^.^ with the superscript 
denoting one of the two types of boundary conditions described. As far as infinite volume limits 
are concerned, there are two alternative scenarios: (1) for fixed JJ take A /" 1? and (somehow) 
average the resultant measure over J. (2) Take the (vague) infinite volume limit of the finite volume 
quenched measures described above. It turns out that for the pertinent cases at hand, the order of 
the procedures is immaterial and the result is independent of how A /* 1? . We will provide a brief 
sketch of how this is established in the appendix. 

Below we summarize our claims concerning these quenched measures. 
Proposition Al Consider the Hamiltonian in Eg. 113. For all F satisfying the ferromagnetic 
condition, that is j > with probability one, the limiting measures 

Q 1 f ] =}%Qf\a (3-3) 

exists (independently of how A /" Z 2 ) and similarly for QjP . Furthermore, the measures and 
QjP are invariant under Z 2 shifts - as well as other Z 2 symmetries. 

Proposition A2 Consider the Hamiltonian in Eg. 13.1]) . If fi^ and /ij' are limiting Gibbs measures 
(which also exist independently of how the infinite volume limit is taken) then the average of these 
measures is well-defined and equal to their respective Q counterpart. Finally, if F is such that 
the spontaneous magnetization as defined thermodynamically ( or as will be discussed subsequent to 
Eg. 13.9]) ) is zero then the limiting quenched measure - satisfying all properties of Proposition Al 
above and A3 below - is unique. 
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Proposition A3 The limiting measure Q F satisfies the strong mixing condition. So, in particular, 
if the magnetization vanishes the unique measure is strongly mixing. 

While most likely these are non-Gibbsian measures, as is evidenced by the result in |13j on 
a related system, they are physically motivated and, possibly, experimentally accessible. Notwith- 
standing, at self-dual points (and presumably at other critical points) the combined efforts of species 
2 through q will fail to achieve percolation in spite of their high density. Furthermore in the dilute 
case (and related cases) despite a low density, the type-1 spins can, on their own, achieve percolation 
in the Q^;' - measures. 

3.2 Statement of Main Results 

For < J < oo, the dual coupling is defined by 

J* (J) = log + (3.4) 

and the dual model is defined by the assignment of the coupling J*(Jij) = J*j to the bond 
of the dual lattice (Z + i) 2 that is traversal to the bond (i, j). A model is self-dual[f| if 

dF(J id ) = dFtfJ. (3.5) 

Our result on the bond-random Potts ferromagnets reads as follows: 

Theorem 3.1 Suppose that for the ferromagnetic bond strength distribution F , in both the direct 
and the dual model, the spontaneous magnetization vanishes. Then, with probability one, in the 
(unique) limiting quenched measure there are infinitely many circuits surrounding the origin such 
that the spin-type is constant throughout the circuit and, moreover, there are infinitely many such 
circuits of each spin-type. In particular, there is no percolation (or even * -percolation) in the 
various marginal measures which identify as many as (q — 1) of the spin-states as a single state 
notwithstanding that the density of this amalgamation is 1 — - . 

Corollary 3.2 For a 2D disordered Potts model at a point of self-duality, the (hypotheses and) 
conclusions of Theorem \3.1\ hold. 

2 For models that do not respect all of the Z 2 symmetries, a more general definition of self-duality is possible 
e.g. the vertical bonds distributed as the dual of the horizontal bonds and vice versa. While most of our results go 
through easily in these cases, we make no further reference to these extensions since the present purpose is to construct 
"natural measures" on Z 2 
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Our next result concerns situations which do have percolation: 

Theorem 3.3 Consider the dilute model with parameters a and f3 defined by J t j = (3 with proba- 
bility a and zero otherwise. If a > h there is a f3 c < oo such that for all f3 > (3 C , the magnetization 
m((3,a), given for [1] boundary conditions by the excess density of species 1 above |, satisfies the 
inequalities 

< m(/3. a) < Poo (a) (3.6) 

where Poo (a) is the percolation probability in the independent bond-model on 1? . In particular, for 
all [3 > (3 C , in the measure Cfp' , there is (with probability one) an infinite cluster of species 1 while 
the overall density of this species is m(J3, a) + -| < Poo (a) + g> by considering large q and a close to 
this density can be made as small as desired. 

Both features may be exhibited in a single model: 

Theorem 3.4 Let a G (0, 1) and K denote a bounded random variable satisfying K > b with 
probability one with b considered "large". Consider a disordered Potts ferromagnet with strong and 
weak bonds: Suppose that for each {i,j), with probability a, Ji j is equal, in law, to K while otherwise, 
with probability (1 — a), it is equal in law to J*(K). Then for b large, the magnetization is very 
nearly Poo (a) and the overall density in the measure Q^p is very nearly Poo (a) + |. Since the model 
is self-dual at a = |, both the high and low density percolative phenomena occur as the parameter a 
passes through \. 

3.3 Graphical Representations, Dual Models and Dominations 

For (i,j) a neighboring pair let us define 

Ri,j = c Jl - 3 - 1. (3.7) 

As is well known |f4j the model admits the random cluster representation: In finite volume, if lo is 
a configuration of occupied and vacant bonds (or edges) the probability of lo is given by 

P» ;A Hcx [] ^ C * M (3-8) 

where (i,j) £ ui represents the event that the particular bond is occupied and cr{ui) denotes the 
number of connected components which are counted by rules (jj) in accord with conditions specified 
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at the boundary: ft = [1] and (j = [f]. In the former case - sometimes called the wired measures - 
cW (cj) counts all clusters that are attached to the boundary as part of the same connected component 
while c[ f l(oj) simply counts the number of components in u> by the conventional definition. Given 
these random cluster measures, it is possible to write down the conditional probabilities of spin 
configurations given a bond configuration. This is done by insisting that the spin-value is constant 
throughout each component of lo and, except for the components attached to the boundary, assigning 
the spin-types to each component independently and with equal probability. As for the boundary 
component, in so far as concerns the two setups in this work, the procedures are simple: For [f], 
nothing special is done - all boundary components are treated like the internal components. For the 
wired or [1] case, all components of the boundary are set to the first spin-state. Thus, for example, 
the finite volume magnetization at the origin for specified J is given by 

m J;A (0) = Eg]v(W) ~\ = p ji1({0 - dA}) (3-9) 

where {0 <-> dA} is the event that the origin is connected to OA by an occupied path. The average 
magnetization at the origin mj^(0) is obtained by averaging (|3.9p over F. The infinite volume 
spontaneous magnetization mW(0) is obtained by taking the limit A f Z d of mj^(0). We call 
mW (0) the spontaneous magnetization of the system. 

To conclude: on the basis of Eq. (|3.9p the magnetization at the origin is non-zero in any given 
quench (realization) if and only if the right hand side does not tend to zero as A f Z 2 ; i.e. percolation 
and spontaneous magnetization are synonymous. That the limiting percolation density exists (with 
or without the quenched average) is well known and anyway can be derived on the basis of what 
is discussed later. The fact that percolation probability is equal to the thermodynamically defined 
spontaneous magnetization has been proved elsewhere, see |16j and references therein, as well as [4] 
and [3J. 

In finite volume, the dual model is defined as follows: If is an edge of Z 2 , let denote 
the corresponding edge of (Z + |) 2 and let 

R*, r = (3.10) 

which is the equivalent to Eq. (|3.4[) via Eq. (|3.7[) . For a finite A C Z 2 - with A regarded as a graph 
- consider the dual graph, A* consisting of all (dual) edges corresponding to the (direct) edges in 
A and the collection of (dual) sites which are the endpoints of these dual edges. As is well known, 
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the model on the dual graph with parameters R**j* has configurations which are in one-to-one 
correspondence with (and have the same probabilities as) the original setup on A. Of course, some 
attention must be paid to the conditions at the boundary. All that is needed in this work is the 
readily verified fact that the model with wired boundary conditions on A associates with the model 
with free boundary conditions on A* and vice versa. We will refer to the initial model as the direct 
model and the induced distribution for the J?„ •» - collectively denoted by I* - by F* . Of course 
when it comes to integration, we may use the "direct" dF. 

A model is said to be self-dual if the probability distribution of the {J*j} is the same as that of 
original, e.g. Eq. (|3.5p . or equivalently, 

Rid =d (3.H) 

The "nicest" examples concern a self-duality which is achieved according to a temperature param- 
eter, in which case one can speak of the self-dual temperature (3, c.f. Eq.(30) in [12] , but this will 
not be necessary in the present work. 

For the sake of completeness, let us recapitulate in brief (special cases of) the domination argu- 
ments that were derived in [H[3]. For a fixed bond (i, j) in finite volume with free or wired boundary 
conditions, let us calculate the conditional probability that the bond is occupied. It is not hard to 
see that this is exactly 

PU = T~~r> (3-12) 



1 + R 

if {i <-> j} while the probability is 



rfS = ^Nr = —ttt \ (3-13) 

q + Ri,j Pi,j + q[l - Pi,j) 



if the endpoints are not connected. (In the former, the number of components is unaffected while in 
the latter, it is reduced by one.) It therefore follows from elementary considerations that for each 
quench, the random cluster measures - wired or free - dominate the independent bond measures 
at parameters pf^ and are dominated by independent bond measures with parameters Pij. The 
ease of calculating quenched averages of independent bond measures (percolation on percolation) is 
what lead to the asymptotically sharp results of [3j. Indeed, as is already seen, when the relevant 
^-parameter is large compared with q, the upper and lower estimates do not differ by much. 
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3.4 Proofs of Main Results 

We start with a preliminary result which is most of what is needed for the proof of Theorem 13.11 

Proposition 3.5 Let F denote a distribution of couplings such that the spontaneous magnetization 
of the dual model is zero. Then with probability one for both Q^p and QjP there are infinitely many 
circuits surrounding the origin such that, within each circuit, the spin-type is constant. 

Proof. Let e > and let V C Z 2 denote a finite set containing the origin. Let A D V and T> VA 
denote the event 

T> VA = {cr | 3 a circuit of constant spin-type separating dV and dA} (3-14) 

where by circuit it is meant a graph theoretical "cycle" of vertices. Then, as we shall see, it is 
sufficient to prove (for arbitrary fixed V and e) that Q F (V VA ) > 1 — e whenever A is sufficiently 
large, and similarly for Q^. We let V - and A - denote sets similar to their tilde counterparts 
enhanced by a layer or two at the boundary to avoid discussion of inconsequential provisos caused 
by discrete lattice effects. Let A* and V* denote the dual sets and let T = \ gy* \ ■ Since the 
quenched magnetization in the dual measure vanishes by hypothesis, for every i* 6 dV* the average 
magnetization at i* in the one-boundary conditions in A* (for the dual model) tends to zero as A* 
gets large. Thus, in a large enough volume, for all such i*, 

m F , ;A ,(i*) = J dFE$. A ,(o-i* —-)<¥. (3.15) 

However, according to Eq. (|3.9p , the integrand is the probability, in the dual model, of a dual con- 
nection between i* and dA* . Thus 



f dFP [ j}. A ,({dV* <-> <9A*}) < f dF 



<\8V*\T<e. (3.16) 



However, 1 — Pj». A »({9V* «-> dA*}) is the probability in the direct model of an occupied circuit 
separating dV from dA in the transformed boundary conditions. So we may write 

J dF P^. A ({3 occupied circuit separating dV and dA}) > 1 — e. (3-17) 

But each spin realization associated with such a random cluster event has a circuit of the stated 
type and we have obtained the desired result in finite volume with free boundary conditions: 

Qf ; a(^a) > 1 - e - ( 3 - 18 ) 
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The result for infinite volume follows by monotonicity of the integrand in Eq. (|3.17p in the system 
size: As discussed in the appendix, for free boundary conditions, the integrand increases when we 
consider the same event (involving dA) in a volume A' D A. For wired boundary conditions, the 
result follows because the in any volume, the integrand increases if we replace [f] by [1]. 

Now, to see that the above implies the presence of infinitely many circuits, let (A^Vfe) denote 
a sequence of shapes as above (which exhaust the lattice) and satisfy Ak-i C V& C . . . such that 
the probability of a circuit in the fc th pair exceeds 1 — tk with J2k 6k < 00 ■ By the Borel-Cantelli 
lemma, with probability one, only a finite number of the prescribed circuits fail to appear. □ 

Proof of Theorem \3.1\ Since the dual magnetization vanishes we already have, according to the 
previous proposition, the part concerning the infinitely many circuits. What is lacking is a proof that 
these circuits are not not all of the same spin-type (as would indeed be the case if the spontaneous 
magnetization were positive). Now, using the fact that the direct magnetization also vanishes, we 
may demonstrate, along the lines of Equations (|3 . through (|3.18[) that outside any finite V but 
inside a sufficiently large A the average probability of observing the dual of a circuit composed of 
vacant bonds is close to one - even with (direct) wired boundary conditions. 

The upshot, when both magnetizations vanish, is that even in finite volume, with high probability 
many circuits of dual bonds and many circuits of direct bonds surround the given V. Since both the 
direct and dual circuits can be "freshly constructed" at an increasing sequence of scales, it follows 
that some portion of these circuits separate each other. We are now, more or less, in the same 
situation as Theorem 12. II for the DaC model: There are many direct rings separated from each other 
(and the boundary) by the dual rings. These direct rings are therefore "at liberty" to take on any 
of the q spin-values. In particular for any s G {1, ...(?} and any finite V, 

Q$ A ({3 a circuit of spin-type s surrounding V}) — ► 1 (3.19) 

as A / Z 2 . The result also holds in Q^! A (which anyway leads to the same limiting measure as 
we show in Proposition A2). The conclusion is that with probability one, in the limiting measure 
there are infinitely many circuits of all types surrounding the origin. This establishes the first claim. 
It also establishes the absence of percolation for any amalgamation-alliance since an infinite path 
(connected or ^-connected) of the alliance without the s th spin-state, starting from the origin, is 
prevented by any of the s th state's circuits, the presence of which has probability one. Finally, as is 
evident from the absence of magnetization, the population density of all species is exactly -. □ 
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Proof of CoroUarv \3.S\ As was demonstrated in Q~2], Theorem 1', the magnetization of the disordered 
2D Potts models vanishes at self dual points. This applies equally well to the dual model. □ 
Proof of Theorem \3.3\ This result is, to the largest extent, contained in [3] so we will be succinct. In a 
given quench, each edge corresponding to a zero coupling bond adds nothing to the ferromagnetism. 
The active bonds must be "reoccupied" (i.e. occupied in the random cluster problem on this media) 
with a conditional probability bounded above by R/[l + R] and below by R/[q + R] where R = 
e 13 — 1. The fraction of these reoccupied bonds that belong to an infinite cluster (in the limiting 
[1]— state) constitute the magnetized fraction. After performing a quenched average, it is seen that 
the magnetization is positive whenever 

W > \ (3 ' 20) 

but is bounded above by foo(-qr^) < Poo{o). This is small and positive for all (3 large enough and a 
close to i. The density of ones, in the one state is just this magnetization plus the ambient -. □ 

Proof of Theorem \3.4\ For simplicity we will treat just the case where with probability a, R = 



e 



1 = e b — 1 = B and with probability (1 — a), R = q/B. Obviously, for a = i, the model is 



self-dual. To prove the remaining statements, we reiterate that 

Poo(p cS ) <m(B,a) <Poo(p) (3.21) 

where, in general, p = J dF[R/(q + R)], etc. Thus it is sufficient to find B and a such that both 
p eS and p are in slight excess of h . In this simple case, we may write 

P eB = a^- q + (l-a )T ^ (3.22) 

and 

p = a + (l_ )_i_. (3.23) 

y 1 + B y ! B + q y ' 

It is readily seen that when a > | and B ^> q the desired conditions are met. □ 



4 Discussion 

We have shown here that for the disordered ferromagnetic g-state Potts model on Z 2 at points 
of self-duality there is no percolation for any amalgamation of q — 1 components in the unique 
translation invariant mixing measure Qf- This transpires despite the fact that the density of these 
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combined q — 1 states is 1 — q^ 1 which will be arbitrary close to 1 as q — > oo. Of course we expect 
that the above holds at all critical points of these models. In particular, we suppose, on general 
grounds, that the conditions of Theorem 13.11 ( namely no percolation and no dual percolation in the 
graphical representation) holds at any critical point in a model of this sort. 

It is further expected that as soon as the temperature is decreased below the critical temperature 
or the overall interaction strength increased, the magnetization will rise continuously from zero. If 
this scenario is correct, i.e. if the magnetization (like the energy density - which in this context may 
be taken to mean the bond density) is continuous for this system, it would follow that in each of 
these states we would have percolation at a small density. 

Unfortunately, even aided with self-duality, we cannot make such an assertion. There are however 
some physical arguments supported by simulations which indicate that this is indeed the case, 
see [22j [9] and references therein. We may use duality to prove the weaker assertion that the 
magnetization is not already positive at a self-dual point (c.f. the discussion in [5]). Indeed, this 
would imply the identical circumstance in the dual model - with appropriately modified boundary 
conditions - and the latter would be a non-percolative state from the perspective of the direct model 
which would mean the existence of two states with differing energy density which is ruled out by the 
result of [3]. More succinctly, there cannot be a point where the model exhibits a percolative and 
a non — percolative state - so the magnetization indeed vanishes at a self-dual point. But this does 
not quite rule out the possibility of an interval of critical behavior - with vanishing magnetization 
- surrounding a self-dual point with a discontinuity in the magnetism at the endpoint and no 
coexisting non-percolative state at the point of discontinuity. 

Trivial examples of discontinuous order parameters coinciding with critical transitions are abun- 
dant in short-range ID systems at zero temperature e.g. the Ising model. The energy is continuous 
as T J, but m(0) = 1. More intricate examples can be found, here we mention two. First there 
is the (reinterpretation of the) mean-field fc-core transition [26j where, in the presence of two dis- 
tinct divergent length scales - and susceptibilities - there is a discontinuity in the order parameter. 
Second, we mention the well known Thouless effect [35] for the ID Ising model with ferromagnetic 
pair-interactions that decay like the inverse square of the separation. Here, as was proved in [1], 
the magnetization is discontinuous at the transition point but this point is critical in the sense of a 
divergent length scale, a divergent susceptibility and, it is presumed, a continuous energy density. 
It is interesting to note that, from the low temperature side, this transition is indeed the endpoint 
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of a critical phase [TUl [H] • 

In this note we have circumvented these mathematical intricacies and demonstrated percolation 
at small densities by considering dilute and "nearly dilute" models. In the former case, the mag- 
netization is always small and in the latter case, it is demonstratively small in the vicinity of the 
critical point. This certainly does not settle the issue of magnetic/percolative continuity in these 
models, but it does, perhaps, bring us a small step closer. 

5 Appendix 

Proof of Proposition Al The first crucial ingredient in what is to follow is monotonicity of the 
associated (fixed coupling) random cluster measures in finite volume. In particular, if Ai D A2 and 
we compare the wired measure on Ai with (the restriction of) the wired measure in A2 - with the 
same Jjj's in the common territory - then the smaller system FKG-dominates (meaning that both 
FKG measures satisfy the stated stochastic domination). For free measures the situation is similar 
only the domination goes the other way. The second ingredient is that expectations of spin-functions 
may be expressed as expectations of random cluster functions which may be easily decomposed into 
increasing and decreasing functions. For example if A\ C Z 2 denotes a finite set and 1a x denotes 
the indicator of the event that all spins on Ai are of type one, it is not hard to see [4] that 

Bf A (lAj = E| ] A(^ (Al) ) (5.1) 

where K(A) is the number of distinct components that intersect A. Observing that K{A\) is a 
decreasing function (and that 1 /q is less than one) the right hand side of Eq. (|5.ip is the expectation 
of an increasing function. Notice then that monotonicity immediately implies the existence of a 
limiting probability for this particular cylinder, independently of how the A's go to T? . 

For multiple sets, the situation is only slightly more complicated. E.g. consider 1 J 4 1 1 J 4 2 (where 
1,4,, is the indicator of the event that all spins in Ai are of type two). Here, the product of the two 
appropriate random cluster functions is almost the right answer except we must insist that there is 
no connection between A\ and A 2 since the latter would demand that some sites in A\ would take 
the same spin-value as sites in Ai. However, this constraint may be written as 1 — 1{a 1 «a 2 } and, 
we are again in a situation where we may analyze (sums and differences of) increasing functions. 
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The general formula for arbitrary cylinder function is seen to be 

<k(n i*) - <a(i- ( ^ k{a > ] nt 1 - HA^ Ak} \) (5.2) 

and we may use monotonicity in volume to demonstrate the existence of a limiting free measure 
independent of how the limit is taken. 

For the [1] boundary conditions, the overall situation is similar with just a couple of modifications 
that need to be made. The first difference is that now the sets A2,...A q must not be allowed a 
connection to the boundary hence the function in Eq. (l5.2|) is be modified by the insertion of the 
product IIj>i[l — l{Aj<->dA}]- The second adjustment concerns the interpretation and counting of the 
number of clusters intersecting the set A\. Clearly any sites in A\ that are connected to the boundary 
are "already" in the 1-st ate and, needless to say are considered part of the same cluster. Thus, 
q-K{A-L) mus t be replaced by q~ KwtyAl ^[i-{A 1 ^dA}" + <7l{Ai<-<-9A}] where Kw denotes the number 
of clusters counted according to the wired rules. Note that the second factor may be rewritten as 
l)l{Ai<->aA}] which, along with q~ Kw ( Al ), is manifestly increasing. Moreover, the (increasing) 
connectivity events {Aj <-> <9A} have the property that if A D A' then {Aj <-> <9A} c {Aj <-> <9A'} 
which is in accord with the decreasing tendency of the measures in increasing volume. Thus, a limit 
emerges again by means of monotonicity properties. 

It is clear that both of the limiting measures are translation invariant - as well as enjoying all 
other Z 2 -symmetries. For example if Ba is a local event which depends only on the spins in A and 
T is some I? translation, it is manefcst that 

Qf!t(A)( B t ( a)) = Q£! a (Sa). (5.3) 
However, by the existence of limiting measure (independent of how A /" J?) the left hand side 

rri m 

converges to QV(£>t(A)) while the right side tends to Q f (Ba) which is the desired invariance. 
Other invariance properties and these properties for the other measure follow from an identical 
argument. The claims of Proposition Al are all established. □ 

Proof of Proposition A2 We start with a proof of the second statement (uniqueness of the quenched 
measure if the magnetization vanishes) since we will use these arguments in the subsequent proofs. 
We consider a sequence of volumes (A) with boundary conditions which may depend (measurably, 
of course) on J. Now with J fixed, the measure obtained by any particular spin-configuration at 
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the boundary will lead to some form of a random cluster measure. The details need not concern us 
here except to say that all such measures (and therefore combinations thereof) are dominated by 
the wired measure. We again consider some local spin-function which we decompose into a sum of 
monotone random cluster functions. We let V\ C V2 C A with V2 large compared to V\ and A much 
larger still. Since the magnetization vanishes, using the argument of Proposition [375] we are assured, 
even for wired boundary conditions on A, that a circuit of dual bonds separates dV\ from dVi with 
probability tending to one as V2 gets large. This is a negative (decreasing) event so the statement 
also holds in the measure coming from the arbitrary boundary condition - which we now permit 
to tend to an infinite volume limit in any fashion. But meanwhile in the vicinity of V± and V2, as 
far as the random cluster problem is concerned, the situation inside the outermost ring separating 
dVi from dV2 is equivalent to free boundary conditions on this ring. Of course the statistics of the 
location of such a ring will depend in detail on J but given such a ring, the conditional measure 
inside Vi is dominated by the one with free boundary conditions on V2 and dominates the one with 
free boundary conditions on V\. The conclusion is that with large probability the random cluster 
expectations lie between the free boundary condition problems on V± and V2 and we may safely take 
quenched averages. By allowing V2 and then V\ to exhaust Z 2 we have proved any infinite volume 
measure is the same as the limiting infinite volume free measure. 

To prove the first statement, let us distinguish two cases depending on the vanishing or non- 
vanishing of the spontaneous magnetization. When there is no magnetization, the argument is only 
a slight reworking of the discussion in the previous paragraph. First, by appealing to stochastic 
(FKG) monotonicity in volume, arguments similar to ones already used show that for a.e. J, the 
limiting /ij is unique, independent of how the limit is taken, etc. (Thus, we are technically proving 
a stronger result than the one claimed.) So, in particular, for a.e. J, the infinite volume measure 
may as well be considered as the limit of measures constructed in finite volume with free boundary 
conditions. Now consider, for fixed J, the thermal average of some local spin event A. As before, 
this may be expressed as a sum involving expectations of increasing random cluster functions. Let 
V\ denote a large volume which contains the support of the spin function and V2 containing V± larger 
still. With the vanishing of the (quenched) magnetization, again using the arguments in the proof of 
Proposition 13. 51 we are assured that with high probability (in the random-coupling/random cluster 
joint measure) there is a dual ring of vacant random cluster bonds separating the boundaries of 
the V's; in particular, with probability tending to one as V2 tends to Z 2 . This accomplishes two 
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tasks: first it establishes J-measurability of these thermal expectations and second it shows that 
their quenched average is close to the finite volume quenched average with various free boundary 
conditions. By allowing the volumes to exhaust 1? we have established interchangeability of the two 
limiting procedures. 

When the magnetization does not vanish, the argument is somewhat more intricate. To start 
off, we again claim that for fixed J, various infinite volume measures exist. In particular, here we 
will be interested in the limiting wired and free random cluster measures, Pj W ' and Pj f ' that emerge 
from sequences of the corresponding finite volume measures. Furthermore, we will have brief need to 
refer to the quenched limiting random cluster measures that arise from these boundary conditions. 

Let Vk denote any standard sequence of volumes tending to Z d , e.g. nested squares centered at 
the origin. Let O/c denote the event that there is a ring of dual bonds separating dVk from infinity. 
Finally, let denote the event that infinitely many of the 0fc-events occur. We claim that if the 
magnetization is positive, then the set {J | Pj W ' (0) ^ 0} has F-probability zero. Indeed, let us first 
establish that this event is is a tail event. Foremost, it is manifestly shift-invariant. Next, supposing 
that J is such that Pj W '(0) ^ 0, let us demonstrate that PjT'(0) ^ if J and J' differ at only a 
finite number of bonds. This follows easily since, if I? is a finite set of bonds and J = J' on D c , then 
for any fixed configuration on D, the conditional measures agree. Therefore, letting D denote the 
set where JJ and J' disagree and letting D$ denote the event that all bonds of D are vacant, we have 

p [ ;\e) > p^iD^i&iD,) = p[: ] p )p[ wl (©|z? ) > p [ ; ] (D,)p^ ] (®) (5.4) 

where the last inequality is by the FKG property. Noting that by Eqs. (|3.12p and (|3.13[) , PjT^-Dg) 
is not zero (even in the dilute case) the claim follows. So now we must rule out the possibility that 
{J | Pj W '(0) 7^ 0} has F'-probability one. Assuming that this is the case, let us demonstrate the 
implication that there is no percolation in the random cluster version of the quenched free measures. 
To this end, let & e denote the event that a separating ring occurs inside Vg. Under our assumption 
we have, for almost every JJ, 

hm Pj, fl (©^|0) > Pj, w] (0 £ |0) = 1 (5.5) 
Letting e £ R + there is therefore a (random) I such that 

Pf(0 f |0)>l-e (5.6) 
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and while this I may have slow decay, it is clear that it is not defective. Thus, for any small 77, we 
may define an Lq such that £ < Lq with ^-probability in excess of 1 — However, it is clear that 
for every J and any L and L' with V > L, 

Pf(& L \&) = Pf(® L \® L .) < Pf AL ,(& L ) (5.7) 

because free boundary conditions enhance any decreasing event inside. Hence, for any sequence of 
volumes tending to infinity with free boundary conditions, the quenched probability that the origin 
is connected to 8Al is less than e + 7/; i.e. there is no percolation in the limiting version of these 
measures. Meanwhile, there is percolation in the wired rendition of this measure since, by hypothesis, 
the magnetization is positive. A priori the two measures are stochastically ordered and so, evidently, 
the ordering is strict. By (the corollary to) Strassen's Theorem, [37] - c.f. the discussion in [S3] on 
page 75 - this implies that the one-dimensional marginals, which here amounts to the translation 
invariant FK-bond density, differ strictly. Back in the spin-system, by an argument to be found 
in [24] - presumably known to others - an exact formula can be presented that relates the bond 
density in the random cluster models to the energy density in the corresponding spin-system. The 
implication is a difference among states, (and hence a discontinuity) in the energy density. This is 
forbidden by the result of [5] where it was shown that for the systems described by the interaction 
in Eq. (|3.1| . the energy density is continuous. Thus, henceforth, we may safely assume that with 
F-probability one, Pj, w] (0) = 0. 

Now consider a local spin function which, as before, may be expressed as a sum of terms involv- 
ing increasing random cluster functions. Notice that in this case, for the modification of Eq. (|5.2[) 
discussed shortly after the display, the condition that certain sets must not be connected to the 
boundary are now replaced by the condition that these sets are not connected to the infinite cluster 
which is here identified as a cluster of spin-type 1. We may assume the event C so that if V\ (which 
contains the support of the spin function) is large, we are assured that dV\ is connected to this in- 
finite cluster with probability close to one. Further, if V2 which contains V\ is sufficiently larger 
still, we may use the absence of dual percolation to generate, with high probability, an additional 
circuit of occupied bonds separating dV\ from dVi . By courtesy of the infinite path emanating from 
dVi, any such ring represents, from the perspective of the spin-system, sites which are a constant 
spin-value of 1. Thus, by consideration of the measure conditioned on the outermost separating 
ring, we have that with high probability, the random cluster functions have expectations which are 
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bounded between expectations in the wired measures on dV\ and the wired measure on dVi. And, 
importantly, in both measures, the boundary conditions can be identified with spin boundary con- 
ditions of the 1-type. The remainder of the proof now follows along similar lines as the previous 
case. □ 

Proof of Proposition A3 Let / and £ denote two local spin functions which are assumed to be of the 
form appearing on the left hand side of Eq. ()5.2() . We let T denote some Z 2 translation operator and, 
for notational simplicity, let g = T(£) denote the translated version of £. The sets corresponding to 
g will be denoted by B\, ...B q and we will use A's for the function /. Since, eventually, |T| — > oo it 
may be safely assumed that these A-sets and the B-sets are all disjoint. 

While the functions / and g individually admit an expansion as on the right hand side of Eq. (l5.2| . 
the product of / and g docs not expand into the product of the expansions. Indeed, what we obtain 
is exactly the form of the expression in Eq. (|5.2p with the sets Aj replaced by Aj U Bj . Because 
of these - and additional upcoming - difficulties, we shall segregate the high and low temperature 
cases (meaning vanishing or non-vanishing of the magnetization) . We start with the simpler case in 
which the magnetization vanishes and thence, according to Proposition A2, we might as well assume 
that have constructed the limiting measure from free boundary conditions. 

Let V\ denote a box that contains all the A-sets and, as previously, V% a much larger box 
containing V\ with dV\ and dVi well-separated. Similarly let V{ and denote the corresponding 
sets (with the similar characteristics) for the function g. It will be assumed that T is large enough so 
that V2 and V 2 ' are disjoint. We let T and Q denote the random cluster expression for the functions 
g and / - i.e. the term inside the expectation on the right side of Eq. (|5.2p - and K, the expression 
for the product fg. Consider the event R, that a circuit of vacant bonds separates dV\ from dVi 
and similarly for the event R! . It is claimed that should both (or either) of these events occur then 
JC indeed factors into TQ, e.g. 

K.l RnR , = [T][g]l RnRl . (5.8) 

Indeed, under this condition, there cannot be a connection between V\ and V{ so we have K(Aj U 
Bj) = K(A J )+K(B : j) and similarly, the event {AjU-B; <-> Aj\JBj} c reduces to {A { <-> Aj} c n{B t <-> 

The reader is reminded that the absence of magnetization implies that with high probability (with 
respect to the joint measure) such rings are likely to occur. If we therefore defer attention from the 
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unlikely event that either of these ring events fail, we are left, for fixed JJ, with the random cluster 
expectation of a sum of the form . \ijCt$j where the a's and (3's are non-negative increasing 
random cluster functions - bounded by 1 - and the A's are of undetermined sign. Let us therefore 
consider, e.g. in free boundary conditions on some much larger A containing all of the V's, the 
generic expectation Ej. A [(o!)(/3)(lflnfl')]- On the one hand, 

Ej ] A [(a)(/3)(lflnfl')] > Ejjjafl - P J;A ({i? n R'} c ) (5-9) 
and, in the above mentioned spirit, let us focus only on the first term. Using FKG, 

E^a/3]>EjJ>]E| A [/?]. (5.10) 
We must now perform the average over J of both sides. We claim that 

/ dF E£ A [a]Egj/3] > j dF E^ A [a] J dF Egj/3]. (5.11) 

rri 

Indeed, the quantity Ej. A [a] is the random cluster expectation of an increasing function which as is 
not hard to check, is an increasing function of JJ and we may again use the FKG inequality because 
the independent measure enjoys this property. 

On the other hand, we may take the expectation on the left side of Eq. (|5.9p and expand according 
to which pair of rings is the outermost pair in V% \ Vl, and, respectively, V 2 ' \ V{. Under this condi- 
tioning, the expectations are independent; both with regards to the fixed J random cluster measures, 
where each of the factors behaves like the free boundary measure at the respective conditioned rings, 
and with regards to the JFs because the regions of interest are disjoint. Moreover, since a and (3 are 
increasing functions, the conditional expectations are bounded above by the expectations with free 
boundary conditions on the boxes V2 and V 2 '. Thus, all in all, 

Ej A [(a)(/3)(W)] < P J;A ({i?n J R'})Ej 1 v , 2 HEj 1 v .,[/3] < E^ [a]E^ v ,\J3] (5.12) 

where the rightmost term stays factored in the quenched average (over J) due to the aforementioned 
independence. In light of Eq. (|5.11[) . Eq. (|5.12p and their immediate consequences, the individual 
pieces of the function have, after expectation, essential upper and lower bounds in terms of the 
product evaluated in infinite volume. The original functions can then be reconstituted and, as A 
and the V's tend to Z 2 (the latter as the magnitude of the translation becomes large) we obtain the 
limiting equality 

km E^ 1 (/TW)-E| 1 (/)E| 1 W = 0. (5.13) 
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The case of the wired/ type- [1] boundary conditions follows a similar tack but is somewhat more 
arduous. The setup is pretty much the same as in the previous paragraphs: |T| is large, V± and V{ 
are large volumes which contain the supports of the two functions and are, in turn, each contained 
in V2 and V% which themselves are disjoint. All of this takes place in wired boundary conditions on 
a volume A that is much larger still; the wired boundary conditions represent spins of type 1. We 
will be considering the expectation of the random cluster function corresponding to the product of 
the spin functions. 

The important simplifying feature is that, with probability close to one, effective boundary 
conditions will again isolate the two volumes. In particular consider the event, 7Z, that a circuit of 
occupied bonds separates V\ from dVi and that this circuit is connected to 9A; similarly we define 
TV . We claim, on the basis of the arguments in the last two paragraphs in the proof of Proposition 
A2 that 

Jpi A (Tlml')dF^l (5.14) 

as first A and then the V's get large. Thus, for all intents and purposes, we may consider the random 
cluster functions under the condition that TZHTV occurs. Noting that these separating rings are part 
of the boundary cluster (and that we employ the low-temperature modifications described shortly 
after Eq. (|5.2p ) we have, for i = 2, . . . , q 

{Ai u Bi <-> dky n (ft n n!) = [{Ai <-» dv 2 y n {b, <-» dv^y] n {K n n!) (5.15) 

where it has been observed that the under the ring condition, the relevant connections to dA occur 
if and only if the connection to the boundary of the corresponding V% or V 2 ' occurs. Furthermore, 
among these sets, all disconnection provisos break down to conditions among the A's and S's alone. 
I.e. 

{Ai ub,« Aj u B y n (n n H') = {A, <-» A 3 y n {B t <-> Bj} c r\(Kn n'), (5.16) 

2 < 2, j < q, i j- Thus as far as the indices 2 through q are concerned, under the condition of the 
ring events, the function is already in product form. Wc finally turn attention to the term involving 
the cluster content of A\ U B\, Here it is claimed that with the "low temperature modifications" 
discussed subsequent to Eq. ()5.2|) and the condition TZdTi', this term reduces to (q- Kw ( Al ) [1 + (q — 
l)lA 1 ^dv 2 ])(<l~ Kv ' < " Bl * > [1 + (q ~ l)lsi^9v 2 ']) i-e this too is in product form. The remainder of the 
proof is similar enough to previous derivations that we need not spell it out in full detail. The key 
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ingredients are the monotone decreasing property of the wired boundary conditions for all increasing 
events and the fact (used earlier in a related context) that the connections to increasing boundaries 
respect this direction of monotonicity. E.g. 

PZv^i - m}) < p!a(Mi <-> dv 2 } \nmz')< Pj^dA, <-> dv 2 }) (5.17) 

where the middle term pertains to conditioning on the "outermost ring" in the region V2 \ Vy ■ □ 

Corollary. For all temperatures, spin functions of the form Ka — YliPA •"•o , i=i an d positive linear 
combinations thereof are positively correlated in the finite volume and infinite volume [1]— quenched 
measures. Furthermore, the 1-state marginals of these measures conditioned on J have positive 
correlations; i.e. these marginal measures have "conditional positive correlations". In particular, the 
above holds in the unique (free) measure associated with the vanishing of the magnetization. 

Proof. For the type of cylinder function mentioned above, in above mentioned boundary conditions 
the associated random cluster functions are actually increasing. Thus, using a double FKG derivation 
as in Eq. (|5.10p - Eq. (|5.1ip . we have the positive correlations. The second statement, while not 
particularly trivial, was proved in (first Lemma) and generalized in |19j . □ 
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